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The motion in the complex plane of the zeros to various zeta functions is 
investigated numerically. First the Hurwitz zeta function is considered and 
an accurate formula for the distribution of its zeros is suggested. Then func- 
tions which are linear combinations of different Hurwitz zeta functions, and 
have a symmetric distribution of their zeros with respect to the critical line, 
are examined. Finally the existence of the hypothetical non-trivial Riemann 
zeros with i?e s 7^ 1/2 is discussed. 
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1 Introduction 



The Riemann zeta function ({s) [0, 0, ^ and its zeros is a source of inspiration in 
the field of quantum chaos f^. It is well known that the imaginary parts of the zeros 
to ( on the critical line Re s = 1/2 have some striking similarities with a spectra of 
an Hermitian operator with time reversal symmetry broken The famous Riemann 
hypothesis states that in fact all non-trivial zeros of ( lie on the critical line and the 
search for the Hermitian operator, and its classical counterpart, has been intensified 
recently [0]. 

With help of powerful computers the Riemann zeros have been explored in great detail 
IP). In this paper we present much more modeste numerical experiments for the zeros of 
Hurwitz zeta function |^ 

C(s, a) = -. ^, < a < 1, a > 1, s = a + it, (1) 

where ({s, 1) = C{s), and for linear combinations of them, 

m 

^(s) = E ciCis^ai). (2) 
1=1 



Important examples of "^{s) are the L-functions P, 0, [TO] for which a; = l/m and 



c/ = x(/)/m'^, are Dirichlet characters modulo m. The motivation for the study is 
that it is often very useful to study how a system behaves as function of a parameter. 



see e.g. chapter 6 in [jTT|. We will, inspired by the similarities mentioned above, use 
the word "spectra" for the imaginary parts of the zeros, not only for ({s, 1) but also for 
C{s,a) and \l/(s). In sect. 2 we consider the spectra for ({s,a) as function of a while 
we in sect. 3 discuss the motion of the zeros for some ^E'(s) with m = 5, and symmetric 
distribution of the zeros. In the latter case the coefficients, q, depend on one or two 
parameters and it is the motion with respect to their variations that is studied. Bombieri 
and Hejhal has shown that in the limit of large t almost all zeros for these symmetric 



functions are simple and lie on the critical line |]T2|. Finally in sect. 4 some results on 



symmetric zeta functions with m > 5 are presented together with some speculations on 
the existence of non-trivial Riemann zeros with Re s ^ 1/2. The numerical calculations 
presented are done with Mathematica on a Sun Ultra machine. 



2 Hurwitz zeta function 

In this section we focus on the motion of the zeros of ({s, a) with respect to a. Since 
along a zero d( = (' ds + / da da = we get the following equation for the motion of 
a zero at s = z 

dz dC(z,a) 1 z-C(z+l,a) 



The second equality in (3) is obtained by derivating (1) with respect to a. For the 
extension of to the whole complex plane by analytical continuation see p. 

Note that when a varies zeros can merge to multiple zeros or go away to infinity but 
never disappear or be created. At multiple zeros C = so the equation of motion must 
be modified, see sect. 3. A spectra for the Hurwitz zeta function can now be obtained 
by regarding t = Im z{a) as an eigenvalue. To get a spectra at a = 1 with unit mean 
spacing one has to use scaled eigenvalues N{t). With N{T) we denote the number of 
zeros in the range < Im z < T, < Rez < 1, and the following asymptotic formula 

T T T 
N{T) = -Log{-) - -, (4) 

ZTT ZTT ZTT 

is used. To obtain a spectra with approximately unit mean spacing in the whole region 
< a < 1 we suggest that (4) should be modified in the following way 

iV(r,„) = |Loj(|)-|-|L<,9(a). (5) 

where the zeros now can be located outside the critical strip. Numerical experiments in 
the range 10~^° < a <1 , < T < 10^ indicate that (5) is very accurate but we have 
not been able to prove it. It is certainly correct for a = 1/2 since 

C(5,l/2) = (2^-l)-C(s,l). (6) 

The zeros of 2* — 1 have a = with spacings At = j^^^ and their number up to 
Im z = T is given by the last term in (5). The formula also seems to be accurate for 
L-functions if m = p, p is a prime number, and a = 1/p is used in (5), see also ||^. 

In fig. 1 a part of the spectra for the Hurwitz zeta function is presented. The 
30:th to the 43:th zero in the range 10~^ < a < 1 are shown. We can observe that 
the spectra share some properties with a generic quantal spectra since there are many 
avoided crossings but also degeneracies. Note that these degeneracies are only in t, 
not in a. However, we see also that the crossings occur in a somewhat peculiar way. 
They take place when N{t{a)) for one of the zeros is strongly decreasing with increasing 
a and then usually more crossings occur for nearlying a-values. This effect is much 
more pronounced higher up in the spectra. It can also be seen from fig. 1 that the 
spectra becomes more equidistant when a — > . This can be understood qualitatively 
since C{s,a) ~ ^ + C{s, 1) for small a-values. The zeros seem thus to move towards 
s = when a — 0. From the analytical continuation of (1) to the whole complex 
plane we know that ({s, 1) has trivial zeros at s = —2, —4, —6, .... From (3) we expect 
that these zeros move to the right on the a- axis when a decreases and eq. (6) gives 
that these zeros have moved two units at a = 0.5. When a —>■ the zeros approach 
5 = 1,0,-2,-4,-6,.... 



3 Symmetric functions ^(s) for m = 5 

Symmetric functions here means functions of the type given by eq. 2 for which the non- 
trivial zeros are located symmetrically around a = 1/2, i.e., if "^{z) = also \I'(1— ^) = 0. 



The symmetric "^{s) which also are L-functions can be expressed as Euler products 
and then it is straightforward to show that 7^ when cr > 1 [Q. The three 

functions of this type with the lowest m- values are ({s, 1), -^{({s, 1/3) — ({s, 2/3)) and 
^(C(s,l/4)-C(3,3/4)). 

For m > 5 opens up the possibility of constructing symmetric "^{s) which are linear 
combinations of L-functions. For these linear combinations there is generally no Euler 
product and zeros can appear outside the critical strip. For m = 5 we can for example 
make the following ansatz 

^5o(s, /?) = ^(C(^, 1/5) - C(^, 4/5) + [3 (C(s, 2/5) - ({s, 3/5))), (7) 

,0 denotes odd since q = —Cm~i in (7). Following the lines at p. 212 in a condition 
for \l/5o(s,/5) being symmetric can be obtained, 

47r ^ 27r - , 27r ^ An. . . 

sin p szn — = p (sm h p sin — ). (8) 

5 5 5 5 

The two odd L-functions with m = 5 have (3 = ±i. In fig. 2 a part of the spectra for 
^5o(s,/?), with P lying on the circle (8), is shown. To obtain this figure the equation 
of motion (3) is used but with a and Cl-^, «) replaced by P and \l/5o(s,/?). In contrast 
to the spectra for the Hurwitz zeta function double zeros now appear frequently and in 
the figure four of them can be seen. It happens when two zeros meet on the critical line 
and move outside a = 1/2 and then again when these two zeros join on the critical line. 
At multiple zeros 'i!'^^{z,P) = so ^ becomes infinite there. To overcome this difficulty 
a small complex constant can be added to the parameter. For each turn on the circle 
(8) in the counter- clockwise sense a zero from the lower half-plane moves to the upper 
half-plane and this is why most zeros in fig. 2 move upwards. It is supposed, but not 
proven, that for the two L-functions at /5 = ±z all non-trivial zeros lie on the critical 
line. It seems hard to gain insight to that problem from the study of this one-parameter 
family of symmetric functions. It becomes more interesting when we now consider to a 
two-parameter family. 

In the same way as above we can try to construct even, i.e. q = Cm-i, and symmetric 
functions with m = 5, 

^5e(s, P) = ^(C(s, 1/5) + Cis, 4/5) + P ids, 2/5) + C{s, 3/5))). (9) 

However, it now turns out that besides the condition, see again p. 212 in 

Att ^ 2tt - 2tt An 

cos— + P COS— = P [COS— + P COS — ) (10) 



P must also fulfill l + /3 = 0so/? = — lis the only possibility, which is a L-function. To 
make a continous change of the spectra, as for \l/5o, we can take a linear combination of 

■^5e with C(S,1), 

Vl/e(s,/5,7) = ^5e(s,/3) + ^C(s,l). (11) 



The symmetry requirement leads to the following equations for the coefficients (3 and 7 



4:77 ^ 27r 7 - , 27r ^ 47r 7. 
cos— + p COS— + - = p {COS— + [3 cos— + -j 
5 5 z 5 5 z 

^ 27r 47r ^ 

l + (5 + l=l{cos— + Pcos— + ^-). (12) 

For (3=1 the first equation is fulfilled for any 7 while the second equation gives 7 = 
1 + ^/h e*^ where 6 is an angle between and 27r. Besides the symmetric functions 
given by (12) there is one more possibility, namely /? = 7 = 1. At this point, where 
cos^ + (3 cos^ + ^ = 0, ^e(s,l,l) reduces to C(s, 1) but it is not possible to continously 
vary this spectra since /3 = 7 = 1 is not a solution to (12). On the circle P — 1,^ — 
1 + ^/5 e'^ we have 

*e(«,^)=(l + ^)C(«,l)- (13) 

When 7 moves counter-clockwise along the circle the zeros of 1 + ^^^-172 (in the following 
denoted as trivial), separated by At = 2%/ Log{5), move upwards with ^ = l/Log{b) 
while the zeros of ^(s, 1) remain fixed. 

Let us now put (3 = 1 + t e''^ and consider when [3 moves on a circle with radius e around 
13 — I'm the counter-clockwise sense. Then eq. 12 gives 

7 = 1 + V5 e^^^ + e^-^e^^ (14) 
< (/? < 27r. With these coefficients we get 

*e(s, e, <^) = (1 + ^ZTri)as, 1) + 9^{s), (15) 

where 

9^{s) = ^-Y^C(s, 1) + as, 2/5) + C(s, 3/5). (16) 

Compared to (13), \E'e(s, e, ip) in (15) have no trivial zeros in the critical strip. For small 
values of e the two cases must, however, be similar, see fig. 3. The difference is that when 
an upsloping " almost trivial" zero meet an " almost Riemann" zero the two zeros bifur- 
cate out in the complex plane or they interact and exchange character on the critical line. 

The equation of motion for a simple zero, z^{(f), to (15) for a fixed positive e is 

dz, ^ d^e{ze,e,^) 1 _ . e'^'^ 1 . 

dip dip ' %{z,,e,^) '^^^--'^ ^5^-1/2 ^^,'^{z,,e,ify ^ ^ 

The fix points of this dynamical system are the Riemann zeros. On the critical line, 
where Re z^ = 1/2, turning points also appear. For large values of e the second term 
in (15) dominates and for the symmetric function 5^5(5) zeros outside the critical line 



appear frequently. If z^o denotes a simple zero of ^'5(5) and if ^z{i~p) — z^{(fi) — z^ the 
equation of motion (17) becomes to lowest order 

i2(p i2ip 

^-'(^) - ^ (^^2 - + 1) (18) 

This linearised system gives rise to straight line motion if Re z 00 = 1/2 and cycles when 
Re Zoo 7^ 1/2. These cycles are counter-clockwise for Re z^o < 1/2 and clockwise for 
Re Zoo > 1/2. Thus, for large values of e the zeros of ^e(s, e, near the zeros of (75(5) 
perform vibrational or circulating motion. However, far outside the critical strip there 
is still an upward flow of almost trivial zeros since ^e(s, e.Lp) ^ 1 + for large values 
of a. When e decreases this upward flow moves towards the critical strip. For e > 2 
the flow from the lower plane to the upper half- plane takes mainly place outside the 
critical strip. Note that ^e("S, e, </?) has a pole at s = 1 except for e = 2, (f = n, which is 
a L-function. This shows that a zero from the upward flow is located at s = 1 for these 
parameter values. For e < 2 the flow from lower to upper half- plane takes mainly place 
inside the critical strip. When e decreases the "train" of zeros moves thus to the left. 
Then zeros of g5{s) with a > 1/2 are approached. If C(-^oo, 1) 7^ these zeros must be 
circumvented. The upward moving zeros can pass by these obstacles by a bifurcation, 
i.e. merge to a double zero, with the zero which form a circuit in the clockwise sense 
around the zero, Zoo, to (75 (s). So for each simple zero 2:00 with cr > 1/2 one extra zero 
can be added to the upward flow. 

There is however another possible obstacle for the moving zeros, namely the famous 
hypothetical zeros, Zq, of ({s, 1) outside the critical line. The linearised motion around 
these zeros is also of the form (18), except for a change of sign, so if i?e > 1/2 the 
circuits are now in the counter-clockwise sense. When the flow of zeros comes in to the 
vicinity of such a zero, Zq, one zero leaves the "train" and starts to circulate around 
zq. Thus, a bifurcation is needed if zq is not a common zero of ({s, 1) and g5{s). The 
parameter values for the bifurcation that "creates" the zero circulating around zq are 
below denoted by and (fB- 



4 Symmetric functions with m > 5 

Symmetric functions of the form (15) can be constructed for higher m values. Here we 

concentrate on primes . For m = 7 there are three complex parameters and with the 
ansatz Pi = 1 + exie*'^, = 1 + 6X26*'^ and 7 = 1 + V^e*^*^ + ex^e'^'^ the symmetry 
conditions lead to the following linear equation for the real parameters X = {xi, X2, xs), 



I 



cos 



COS 
7 



47r 
7 



2 



COS 



COS 



67r 
7 

2-K 



COS ^ 



47r 

COS Y- 



2 



(19) 



The rank, r, of this matrix is two and the null space is one dimensional so there is only 
one possible construction, namely 

*e(s, e, = (1 + T^)as, 1) + e—97{s, x), (20) 

where 

gr{s, x) = Xi(C(s, 2/7) + C(s, 5/7)) + X2(C(s, 3/7) + C(s, 4/7)) + x^Cis, 1). (21) 

Numerical investigations of the corresponding matrices for all prime numbers p < 1000 
shows that p — Ar ±1 and the order of the matrices is so the dimension of the null 
space increases with increasing p. 

An interesting case is p = 13. Here the null space is three dimensional so for each point 
X on the sphere S"^ there is a symmetric function 

*e(s, e, ^, X) = {1 + Y^^)C(s, 1) + e—gi3{s, X). (22) 

Let us as in sect 3 assume the existence of simple Riemann zeros, Zq, outside the critical 
line. To each point X there is one point (e^, (ps) where the bifurcation near zq ( see 
sect. 3) takes place and for —X the corresponding point is {eb^^b + tt). If now X 
moves one circiut on the equator of 5*^ we expect that a closed curve (eB(X), ^^^(X)) 
will be traced out. Let us assume that £5 = 6 does not lie on the curve. If then the 
equator is continously deformed to the north pole , say, the closed curve must for some 
X cross £5 = which is impossible for simple zeros! The situation is different for the 
bifurcations near z^q since the location of z^q in the complex plane depends on X and 
Zoo is generally located on the critical line for some X values. Further work on these 
interesting symmetric functions is in progress. 
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Figure 1: The imaginary part of the 30:th to the 43:th zero of ({s, a) in the range 
10~^ < a < 1 is shown. Scaled values ,N(t, a) of t = Im z have been used. For this, and 
all other figures, a 4:th order Runge-Kutta method has been used to solve the equation 
of motion. The suggested formula (5) for N{T, a) is seen to work well. The somewhat 
downsloping overall structure is probably due to higher order terms in (5). 
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Figure 2: A part of the spectra for the symmetric function ^'5o(s,/3) in (7). The pa- 
rameter j3 Ucs on the circle (8) and an angle is used in the figure to parametrise this 
circle. For (3 = +i and —i, which corresponds to left (and right) edge and (f) ^ 5.73 
respectively, arc L-fuctions and all zeros for them are supposed to lie on the critical 
line. Note the two upsloping segments which join the double zeros on the critical line. 
To remedy the singular equation of motion at the bifurcation points a small complex 
constant has been added to the parameter cf). 



Figure 3: A low lying part of the spectra for \f el-s, e, y?), given by (15), is shown for 
e = 0.01 and < < 27r. At e = the spectra consists of upsloping trivial zeros 
crossing horisontal Riemann zeros. Here these crossings are replaced by bifurcations out 
in the complex plane or interactions on the critical line. As for fig. 2 a small complex 
constant has been added to </? to overcome the bifurcations. 
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